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The set of kinetic equations describing the process of conversion of a beam of protons into mesons 
and then to neutrinos is solved. The asymptotic evolution of the density profile of neutrinos is 
essentially the same as that obtained in the previous version of the note for a simple model of 
uniformly damped wave-packet. It shows again that the recently reported "superluminal neutrinos" 
[l[ could be considered as a purely statistical effect due to the fact that the detected neutrinos 
represent a biased sample of initial protons. 

Consider a one-dimensional model which describes propagation of a beam of particles with the speed of light. The 
beam consists of protons with the density p(x,t), mesons with the density m(x, t) and neutrinos with the density 
n(x,t). The kinetic equations read 
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where r p ,r m ,r„ are damping rates for protons, mesons and neutrinos respectively, while 7 P and 7 m describe trans- 
formation rated of protons into mesons and mesons into neutrinos. 

The system of equations ((T|) can be easily solved and for the initial conditions : p(x,0) = f{x) , m(x,0) — and 
n(x, 0) = yields 
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Under the physically justified condition T n « Tp,r m and for t » l/T n the expression d2j) simplifies to 

n{x,t)^^e-^f{x-ct). (3) 

which essentially coincides with the result of the previous version of this note. For the initial Gaussian density profile 
(concentrated at the origin for t = 0) the neutrino density evolves asymptotically as 

n(x,t) =iVexp{- — [a;-ct] 2 }exp{-rV}- (4) 

where d is the width of the density profile and N - the normalization constant. The time resolved density profile (|4]) 
at the point x = L » c/T n can be recast into the form 

n(i,t) = iVexp{-^(^) 2 }exp{-^[(L-^)-rf] 2 } (5) 



where the shift 
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can be misinterpreted as a consequence of a higher velocity v ~ c(l + 8L/L). The origin of this shift is purely 
statistical, the neutrinos created at the longer distance from the proton source, and hence arriving earlier at the 
detector, are slightly less damped than those created at the shorter distance. 

To apply this result to the OPERA experiment [lj one should notice that in the experiment one compares the 
density profile for the protons - the "grandparents of the neutrinos" with the profile for neutrinos. The former is 



2 



measured under conditions where damping can be neglected and the time-resolved shape does not display the shift 
8L. On the other hand, the highly damped density profile of neutrinos is shifted. The damping rate T n describes 
an overall effect of all processes which reduce the fraction of neutrinos produced by the source in CERN which are 
detected in Gran Sasso Lab. Finally, one can compare the relevant numbers. The initial proton density profile is 
composed of five peaks which can be approximated by the Gaussians of the width d ~ 600?ti. This value combined 
with the value of L ~ 700fcm and the measured relative shift SL/L ~ 2.5 x 10~ 5 yield the overall damping factor (the 
prefactor 7p7 m /r p r m is neglected) 



This value is consistent (within a reasonable accuracy) with the ratio N v /N v ~ 10 16 where N v ~ 10 4 is the total 
number of neutrinos detected in OPERA and N p ~ 10 20 is the number of the corresponding protons. 
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